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We measure the angular dependence of the resonant dipole-dipole interaction between two individual Rydberg 
atoms with controlled relative positions. By applying a combination of static electric and magnetic fields on the 
atoms, we demonstrate the possibility to isolate a single interaction channel at a Forster resonance, that shows 
a well-defined angular dependence. We first identify spectroscopically the Forster resonance of choice and we 
then perform a direct measurement of the interaction strength between the two atoms as a function of the angle 
between the internuclear axis and the quantization axis. Our results show good agreement with the angular 
dependence oc (1 — 3 cos 2 9) expected for this resonance. Understanding in detail the angular dependence 
of resonant interactions is important in view of using Forster resonances for quantum state engineering with 
Rydberg atoms. 
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Among the systems currently considered for quantum engi¬ 
neering with applications to quantum simulation [1] and quan¬ 
tum information [2], ensembles of individual Rydberg atoms 
are promising [3-5], since they provide large interactions [6] 
and are scalable [7]. For that purpose, one needs to control 
pairwise interactions in the system. Significant achievements 
using small numbers of atoms have been obtained based on 
the phenomenon of Rydberg blockade [8-15]. The energy 
transfer between atoms [16, 17] observed in the presence of 
resonant interactions also leads to interesting many-body dy¬ 
namics governed by spin-exchange Hamiltonians [18, 19]. In 
those examples, one has to pay attention to the angular depen¬ 
dence of the interaction [20-23]. 

The anisotropy of the dipole-dipole interaction plays a cru¬ 
cial role in solid-state nuclear magnetic resonance [24], and 
in systems such as cold polar molecules [25, 26] or dipo¬ 
lar quantum gases [27], Although experimental evidence of 
the anisotropy of the resonant dipole-dipole interaction be¬ 
tween Rydberg atoms has been observed in confined geome¬ 
tries [28], to date, a direct measurement of this angular de¬ 
pendence has not been obtained due to the complexity intro¬ 
duced by multiple interaction channels and inhomogeneities 
in disordered ensembles [20, 29]. Here we use a combination 
of electric and magnetic fields to isolate a single interaction 
channel, by tuning two atoms separated by a controlled dis¬ 
tance R to a Forster resonance. We measure the evolution of 
the system [30] to extract the interaction strength as a function 
of the angle 9 between the interatomic axis and the quantiza¬ 
tion axis z. Our results further extend the possibilities offered 
by Rydberg atoms for quantum state engineering. 

The electric dipole-dipole interaction between two atoms 
(k = 1, 2) is described by the operator 
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where dk.± = T(dk,x ± idk, y )/V2 are the components of 
the dipole operator in the spherical basis. The operator T4d 
in Eq. (1) contains terms with angular prefactors A\{9) = 
(1 — 3cos 2 #)/2, A2(9) = 3sin0cos0/-\/2 and A 3 (0) = 
3 sin 2 0/2, that couple pair states where the total magnetic 
quantum number M = m\ + m 2 changes by AM = 0, 
AM = ±1 and AM = ±2, respectively. In the absence 
of external fields, two atoms prepared in the same state gener¬ 
ally interact in the van der Waals regime, where several states 
contribute to second order to the interaction [7]. Resonant 
dipole-dipole interactions occur when two pair states coupled 
by Vdd are degenerate, giving rise to stronger interaction en¬ 
ergies Em oc 1 /R 3 [6, 7]. 

Resonant interactions between Rydberg atoms can be ob¬ 
served by applying an electric field, to reach a Forster reso¬ 
nance [31-38]. Here we study the resonance 

59D3/2 + 59D3/2 <—> 61 P 1/2 + 57P5/2 ( 2 ) 

for two Rubidium atoms prepared in \d) = \59D 3 / 2 ,m.j = 
3/2) [30], When 9 = 0, the angular prefactors A 2 and A 3 
vanish, and we only need to consider dipolar couplings with 
AM = 0. In this case, Eq. (2) describes the transition of 
one atom to the magnetic substate | p) = |61 Pi/ 2 ,mj = 
1/2) while the other atom evolves to the magnetic substate 
|/i) = \57F 5 /2,rrij = 5/2). We define the symmetric state 
b/i)s = (\pfi) + \fip}) /V2. In a zero electric field, the en¬ 
ergy splitting between \dd) and \pfi) s is only h x 8.5 MHz, 
and vanishes for a small electric field (F les ~ 33 mV/cm) 
due to the differential Stark effect between the two states [39]. 
At resonance, the eigenstates of the system are the com¬ 
binations |±i) = (| dd) ± |p/i) s ) /v/2, that have energies 
E± = ± \/20 3 /R :i . The coefficient C 3 , defined such that 
the matrix element of t/jd between | dd) and \pfi) for 9 = 0 
is equal to C 3 /R 3 , was measured previously in [30] to be 
C 3 /h = 2.39 ± 0.03 GHz.^m 3 . 

Here, we focus on the angular dependence of this inter¬ 
action. When 0 / 0, Vdd couples | dd) to different mag¬ 
netic substates \pfi) B , as detailed later (see Fig. 4). Since the 
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states \pfi) s have different polarizabilities, it is possible, by a 
proper choice of parameters, to isolate one interaction chan¬ 
nel [37, 38] showing a well-defined angular dependence. This 
is what we demonstrate here for the resonance between \pfi) s 
and | dd), with full control over the geometry of the system, 
which allows us to measure an interaction strength varying as 

Am- 

Our experiment starts by loading two 87 Rb atoms from a 
magneto-optical trap into two microscopic dipole traps, cre¬ 
ated by focusing a laser beam of wavelength 852 nm using 
a high-numerical aperture lens [40]. We control R and 9 
using a spatial light modulator to create arbitrary trap pat¬ 
terns [41], Using a set of electrodes, we control the elec¬ 
tric field F, at the position of the atoms. We apply a 3.3 G 
magnetic field along z, and we optically pump the atoms in 
| g) = |5<Si/2,-F = 2 ,niF = 2). For the Rydberg excita¬ 
tion from | g ) to |d) = \ 59 D 3 / 2 ,rrij = 3/2), we perform a 
two-photon transition using a 7r-polarized laser beam of wave¬ 
length 795 nm and a er + -polarized laser beam of wavelength 
474 nm. During that time, the microtraps are switched off. At 
the end of the sequence, we turn on the trapping light again 
and shine resonant light of wavelength 780 nm on the atoms. 
We read out the state of the system by detecting the fluores¬ 
cence of ground-state atoms (atoms that are in a Rydberg state 
at the end of the excitation sequence are not recaptured in the 
microtraps). We repeat the sequence ~ 100 times to recon¬ 
struct the two-atom state populations P gg , P gr , P rg and P rr , 
where r stands for any Rydberg state and g stands for the 
ground state. 

In a first experiment, we identify the various resonances 
by measuring Rydberg excitation spectra as a function of the 
electric field for two atoms with H = 9.1 /jin and 9 = 32°. 
For electric fields that vary between F z = 4 mV/cm and 
F z = 44 mV/cm, we shine on the atoms an excitation laser 
pulse of Rabi frequency Q/2tt ~ 0.76 MHz and duration 
t = 1 /is. We record the probability P rr to excite both atoms 
to a Rydberg state as we scan the frequency detuning S of the 
excitation laser with respect to the Rydberg line measured in 
a zero electric field. Figure 1(a) shows the result of this mea¬ 
surement. For electric fields smaller than a few mV/cm, we 
observe a single excitation line showing the transfer of popu¬ 
lation from | gg) to | dd) via the pair states \gd) and \dg). At 
resonance between | dd) and \pfi) s , the state | gg) is partially 
coupled to the eigenstates |±j) = (| dd) ± \pfi) s ), and we 
observe the presence of two excitation lines. Because \pfi) s 
have different polarizabilities, resonances between dd) and 
\pfi) s occur at distinct F z , and we observe different avoided 
crossings between | dd) and \pfi) s . The measured spectrum 
is in good agreement with the calculated spectrum shown in 
Fig. 1(b), obtained by simulating the experimental sequence 
for the system, as detailed later in this paper (see Fig. 4). For 
9 = 0, the same measurement (not shown here) gives a single 
avoided crossing between | dd) and |p/i) s [30]. The observa¬ 
tion of multiple resonances when 9 ^ 0 is a first indication of 
the anisotropy of the interaction. 

We now focus on the resonance at F z ~ 33 mV/cm (red 
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FIG. 1. (color online) Rydberg excitation spectrum for two atoms 
with R = 9.1 gm and 9 = 32°. The double excitation probabil¬ 
ity P rr is represented as a function of the detuning 5 and the elec¬ 
tric field F z . The red vertical dashed line shows the measured posi¬ 
tion of the Forster resonance detailed in this paper. The bold arrows 
show the calculated positions of the different resonances, (a) Exper¬ 
imental spectrum, (b) Theoretical spectrum obtained by solving the 
Schrodinger equation for the system (see text). 


dashed line on Fig. 1) and study its angular dependence. For 
a given angle 9, we obtain the interaction strength between 
the two atoms by measuring the time evolution of the system 
in the presence of resonant interactions [30]. In this experi¬ 
ment, R = 9.1 /tm, fl/27r ~ 5 MHz and 9 varies between 
—90° and 90°. We first place the system out of resonance 
(F z = 64 mV/cm), where the atoms show weak van der 
Waals interactions < 1 MHz. We apply an optical 7r-pulse to 
prepare the two-atom system in the state | dd). We then switch 
on abruptly strong resonant interactions by pulsing the elec¬ 
tric field to 33 mV/cm with rise and fall times < 10 ns, thus 
inducing back and forth oscillations between | dd) and |p/i) s . 
After an interaction time 1, the probability for both atoms to 
be in the state | d) reads 


Pdd(t ) = cos 2 


(y/iCM0)\ 
^ hdp J 


(3) 


where we neglect the influence of other interaction channels. 
At the end of the sequence, we freeze the dynamics by switch¬ 
ing the electric field back to F z = 64 mV/cm. To read out 
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FIG. 2. (color online) Coherent oscillations at resonance for angles 
varying between 9 = 0 and 9 = 92°, where R = 9.1 fin l. The 
typical error on P gg is ±0.05. Lines are fits to the data by damped 
sines (see text). 


the state of the system, we shine a second optical 7 r-pulse on 
the atoms, that couples \dd) back to | gg) while leaving |p/i) s 
unchanged. We measure the probability P gg to detect both 
atoms in the ground state. Assuming perfect excitation and 
detection, P gg coincides with P,id ■ 

Figure 2 shows the coherent oscillations observed for dif¬ 
ferent angles 9. The anisotropy of the interaction is evident on 
the data. For 9 = 0 we observe contrasted oscillations with 
frequency v ~ 8.6 ± 0.4 MFlz, compatible with the value 
of C 3 measured using the same experimental sequence and 
varying R [30]. As we increase 9, we initially observe a de¬ 
crease of v. For 9 ~ 57°, the observed dynamics are essen¬ 
tially described by an exponential damping. This is compat¬ 
ible with the expected behavior close to the “magic angle” 
( 9 m ~ 54.7°) where the strength of resonant interactions van¬ 
ishes (Ai(9 m ) = 0), and where one thus expects Pdd to stay 
constant. Further increasing 9 leads again to the observation 
of contrasted oscillations with increasing frequencies. The ob¬ 
served damping originates mainly from dephasing, due to the 
finite temperature of the atoms and fluctuations in the volt¬ 
ages applied to the electrodes leading, respectively, to shot- 
to-shot fluctuations of R and of the detuning between | dd) 
and |p/i) s . The observation of mainly one frequency in the 
oscillations confirms that the interaction channel of choice is 
well-isolated. 

We obtain the interaction strength as a function of 9 by fit¬ 
ting the measured coherent oscillations with damped sines of 
frequency v. Figure 3(a) shows a plot of // as a function of 9. 
Our measurement procedure does not give access to the sign 
of the interaction, and we thus infer it from the expected func¬ 
tional form. The error bars represent statistical errors from the 
fit. For 9 ~ 57°, close to the magic angle, we do not observe 
any significant oscillation within the limits of experimental 
accuracy. In this case, we assign the value 0 MHz to the fre¬ 
quency, with an error bar representing the upper frequency 
one could possibly infer from the measurement and equal to 
l/(2i max ), where f max = 700 ns is our largest interaction 
time. The data show good agreement with the expected angu¬ 




FIG. 3. (color online) Angular dependence of the interac¬ 
tion. (a) Graph of v as a function of 9. The solid line plots 
—C3Ai(9)/(hR 3 ). Error bars show statistical errors in the deter¬ 
mination of v. The shaded area corresponds to a systematic 5% error 
on the calibration of R. (b) Representation in polar coordinates. By 
symmetry, the points at angles at 9 + 180° are taken identical to the 
points at 9. 


lar dependence oc A\(9) plotted as a solid line in Fig. 3(a), 
where we used the coefficient C 3 measured in [30], The polar 
representation of the fitted frequencies as a function of 9 ex¬ 
hibits the characteristic shape of the function (1 — 3 cos 2 9 ), 
as can be observed in Fig. 3(b). 

In the remaining of this paper, we analyze theoretically our 
results. We first consider the situation where no magnetic field 
is applied on the atoms. Coming back to Eq. (1), terms with 
the angular prefactor A\{9) couple | dd) and |p/i) s (AM = 
0). Terms with the angular prefactor „4 2 ($) couple | dd) and 
I Ph) s , with |/ 2 ) = |57F 5/2 , rrij = 3/2) (AM = -1). Fi¬ 
nally, terms with the angular prefactor A 3 (9) couple | dd) and 
|p/ 3 >., with I/3) = |57 F 5/2 , mj = 1/2) (AM = —2). In the 
absence of magnetic field, we thus expect three resonances 
between | dd) and the states \pfi), lp/2) and lp/3). 

We now take into account the effect of the 3.3 G magnetic 
field applied on the atoms during our experiments, leading to 
energy shifts and to mixing of the different states [37, 38]. 
Figure 4 shows the calculated Stark map for the relevant pair 
states, where the dipole-dipole interaction has not been in¬ 
cluded. It was obtained by calculating the (one-atom) level 
shifts of the states |p), |d), and |/j) from the diagonaliza- 
tion of the total Hamiltonian, including the Stark and Zeeman 
effects [42], The electric dipole matrix elements were ob- 





























4 



F z (mV/cm) 


FIG. 4. (color online) Stark map of the relevant pair states in the pres¬ 
ence of a 3.3 G magnetic field. The colored lines plot the eigenen- 
ergies of the two-atom system as a function of the electric field. The 
color scale represents the overlap between the calculated eigenstates 
|$) and the states|p/i) B , \pf 2 )s, \pf 3 )s- The black dashed line shows 
the energy of \dd). The solid arrows show the positions of the differ¬ 
ent resonances. The red vertical dashed line indicates the position of 
the resonance detailed in this paper. 


tained numerically using quantum defect theory [20, 43, 44]. 
We observe that the magnetic field shifts the pair states ener¬ 
gies due to the Zeeman effect, and mixes \pfi) s to \pfi + f) s 
(i = 1,2,3), where |/ 4 ) = \57F r/a , mj = 5/2), |/ 5 ) = 
\57F r/ 2, mj = 3/2) and |/ 6 ) = |57F 7/2 ,m, = 1/2) (see 
solid lines paired by colors in Fig. 4). Because all pair states 
have different polarizabilities, the system exhibits a discrete 
set of Forster resonances (see bold arrows), at the crossings 
between the solid colored lines and the black dashed line. 
The color bar represents the degree of admixture of the dif¬ 
ferent states and shows that the presence of a magnetic field 
increases the number of expected resonances. 

Based on this Stark map, we finally simulated numeri¬ 
cally the experimental spectrum of the two-atom system of 
Fig. 1(a). For that purpose, we first truncated the two-atom 
basis, retaining only the two-atom states differing in energy 
by no more than 40 MHz from \dd), and then calculated the 
restriction of Vjd to this truncated basis. The result of the cal¬ 
culation, shown in Fig. 1(b), is in good agreement with the ex¬ 
periment regarding the positions and the strengths of the dif¬ 
ferent resonances. For our experimental parameters, we note 
that we resolve three avoided crossings. We focused in this 
paper on the crossing between the states | dd) and \pfi) s oc¬ 
curring at F z ~ 33 mV/cm. In this case, the states \pf 2 ) s 
and \pfs)s are off-resonant by more than 10 MHz, which en¬ 
sures that this resonance is well isolated from other interaction 
channels, as observed on the interaction dynamics. On the 
contrary, in the absence of magnetic field, the energy split¬ 
ting between the different states (< 10 MHz) would not be 
sufficient to isolate the different channels. Using a finite mag¬ 
netic field was thus necessary to isolate a single resonance. 
By reproducing the same procedure with different interaction 
channels, one could in principle engineer other angular depen¬ 


dencies ex A 2 ( 0 ) or cx . 43 ( 0 ). 

The control over the anisotropy of the dipole-dipole inter¬ 
action demonstrated in this work will be useful for the imple¬ 
mentation of quantum information protocols [45] or for the 
quantum simulation of many-body spin systems in large 2D 
arrays of dipole traps [41], It is also a motivation to revisit, 
in the presence of anisotropic interactions, effects discussed 
previously such as the prediction of the breakdown of the Ry¬ 
dberg blockade in the presence of nearly resonant interactions 
[46], the existence of few-body Forster resonances [47] or the 
diffusion of spin excitations in complex ensembles of spins 
with long-range interactions [48, 49]. 
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